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ABSTRACT 

Context. We study large trans-neptunian objects (TNOs) using stellar occultations. 

Aims. We derive precise astrometric predictions for stellar occultations by Eris, Haumea, Ixion, Makemake, Orcus, 
Quaoar, Sedna, Varuna, 2002 TX300, and 2003 AZg4 for 2011-2015. We construct local astrometric catalogs of stars 
complete to magnitudes as faint as R = 18-19 in the UCAC2 (Second US Naval Observatory CCD Astrograph Catalog) 
frame covering the sky path of these objects. 

Methods. During 2007-2009, we carried out an observational program at the ES02p2/WFI (2.2 m Max-Planck ESO 
telescope with the Wide Field Imager) instrument. The observations covered the sky path of the selected targets from 
2008 to 2015. We performed the astrometry of 316 GB images using the Platform for Reduction of Astronomical Images 
Automatically (PRAIA). With the help of field distortion patterns derived for the WFI mosaic of CCDs, we reduced 
the overlapping mosaics of CCDs. 

Results. We derive positions in the UCAC2 frame with 40 mas precision for stars up to the catalog magnitude com- 
pleteness limit (about R = 19). New stellar proper motions are also determined with 2MASS (Two Micron All Sky 
Survey) and the USNO B1.0 (United States Naval Observatory B 1.0) catalog positions as a first epoch. Astrometric 
catalogs with proper motions were produced for each TNO, containing more than 5.35 million stars covering the sky 
paths with 30' width in declination. The magnitude completeness is about R = 19 with a limit of about R = 21. We 
predicted 2717 stellar occultation candidates for all targets. Ephemeris offsets with about from 50 mas to 100 mas 
precision were applied to each TNO orbit to improve the predictions. They were obtained during 2007-2010 from a 
parallel observational campaign carried out with from 0.6 m to 2.2 m in size telescopes. 

Conclusions. This extends our previous work for the Pluto system to large TNOs, using the same observational and 
astrometric procedures. The obtained astrometric catalogs are useful for follow-up programs at small to large telescopes 
used to improve the candidate star positions and TNO ephemeris. They also furnish valuable photometric information 
for the field stars. For each TNO, updates on the ephemeris offsets and candidate star positions (geometric conditions 
of predictions and finding charts) are made available in the web" by the group. 

Key words. Astrometry - Occultations - Kuiper belt: general - 



Send offprint requests to: M. Assafin 

* Tables of predictions for stellar occultations by Eris, 
Haumea, Ixion, Makemake, Orcus, Quaoar, Sedna, Varuna, 
2002 TX300, and 2003 AZ 84 for 2011-2015 and respec- 
tive Catalog of star positions for 2008-2015 sky path 
are only available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http:/ /cdsweb.u-strasbg.fr /cgi-bin/qcat? J /A+ A /515 / A32 1 

** Observations made through the ESO run 079.A-9202(A), 
075.C-0154, 077.C-0283 and 079.C-0345. 

*** Also based on observations made at Laboratorio Nacional 
de Astroflsica (LNA), Itajuba-MG, Brazil. 

' Associate researcher at Observatoire de Paris/IMCCE, 77 
Avenue Denfert Rochereau 75014 Paris, France 



1. Introduction 

The trans-neptunian objects (TNOs) are located beyond 
Neptune in a region where only small body differenti- 
ation is expected with regard to temperature changes. 
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This makes them 4.5 billion year old interplanetary 
fossils of the early stages of the formation of the 
outer Solar System. Recent models of planetary migra- 
tion such as the Nice model (jLevison S z Morbidelli 2003 
IGomes et al. 20041 IGomes et al. 2005]) also indicate that 
TNOs are a sensitive laboratory for the study of orbital 
dynamics. Refining the ephemeris and establishing the fun- 
damental properties of TNOs are thus essential for our un- 
derstanding of the origin and evolution of the Solar System. 
This in turn contributes to the understanding of exoplan- 
etary systems and helps astrochemical and astrobiological 
studies (jBarucci 2008j) . 

The analysis of light curves obtained for a stellar 
occultation observed from many sites allows the shapes 
and sizes of TNOs to be measured with a precision of 
a few kilometers. We can probe their atmospheres or 
place limits on their existence down to the nanobar level. 
In addition, occultations can eventually lead to the ob- 
servation or discovery of close binary companions, satel- 
lites, or debris around the central body. For size deter- 
minations, occultations are far more suitable than indi- 
rect estimations such as those coming from albedo as- 
sumptions or from the modelling of optical, infrared (IR), 
and sub-millimetric observations (fCruiksha nk et al. 2005 
lAltenhoff et al. 2004|Stansberry et al. 2008[ ). From these 
direct size measurements, one can derive more accurate 
albedos and place tighter constraints on the surface com- 
position of the TNOs. If the mass of the body can be es- 
timated from the orbits of detected satellites or by other 
indirect means, the density can be estimated more accu- 
rately and thus the internal composition and structure of 
the body can be inferred far more reliably. The detection of 
atmospheres around TNOs will help to improve our current 
understanding of their dynamics and relationship with the 
surface. Detection of Chiron-like jets ([Elliot et al. 1995)) is 
also a possibility. A more robust characterization of com- 
panions/satellite orbits and rotational periods from stel- 
lar occultations and ground-based adaptive optics observa- 
tions could also improve our models for binary formation 
and collision (|Goldreich et al. 20021 |Weidenschilling 2002 
|Sicardy et al. 2011a| ). 

It is more advisable to select large TNOs as targets 
for stellar occultation campaigns. They usually have large 
apparent sizes (30 mas or more) and thus have a greater 
chance of a positive detection. In addition, theory indi- 
cates that larger TNOs are most likely to possess atmo- 
spheres (Elli ot fc Kern 2003)) . In the case of atmospheres, 
even negative occultation chords observed nearby locations 
with positive records are useful, giving valuable upper lim- 
its to the atmosphere size and profile. Larger bodies also 
have smaller relative size errors. The accumulation of pos- 
itive and negative detections also result in significant im- 
provements in the TNO ephemeris, thus refining the orbit 
and substantially increasing the accuracy of the subsequent 
occultation predictions. Improving the chances of success 
minimizes the money and time spent in the complex inter- 
national campaigns, which usually involves large telescopes 
in well-equipped observatories, smaller mobile instruments, 
and specialized personnel. 

Even as one of the most important and larger mem- 
bers of the TNO population, Pluto has only been probed 
in terms of occultations, prior to 2000, almost exclu- 
sively by the 1988 mutual events with its larger satel- 
lite Charon. The first observed stellar occultation by 



Pluto took place only in 1985, when its atmosphere was 
first discovered (|Brosch 1995)) . The atmosphere could only 
be more extensively observed in the next occultation 
in 1988 (jMillis et al. 1993|) . In 2002, a new stellar oc- 
cultation by Pluto was recorded, but with higher time 
resolution flSicardy et al. 2003| Elliot et al.^003| . reveal- 
ing that there had been a dramatic expansion in the 
atmosphere pressure. Another high signal-to-noise ratio 
(S/N) light cur ve was obtained i n the June 2006 stel- 
lar occultation d Young etP al. 2008), which uncovered ev- 
idence of the stabilization of the atmosphere pressure 
in Pluto for the period 2002-2008. This result is sup- 
ported by seven other unpublished occultations observed 
by our group members in Australia, Brazil, Namibia, 
and Chile with 0.5 m to 4 m size telescopes in 2006- 
2010 (Sicardy et al., in prep.). Regarding Charon, after 
its first recorded stellar occultation ( Walker 19 80), four 
others have so far been observed, one in 11 July 2005 
dSicardy et al 2006) iGulbis et al. 2006|) . one in 22 June 
2008 ( |Sicardy et al. 201 la| ), and two other unpublished 
multi-chord events observed in Brazil and Chile (4 June 
2011) and in Hawaii (23 June 2011). These last two events, 
as in 22 June 2008, were double events in which Charon 
and Pluto occulted the same star. In the 4 June 2011 event 
in particular, both occultations could be recorded at each 
site. 

The population of TNOs in general, have been far more 
poorly studied. Up to the time of writing, only seven stel- 
lar occultations have been successfully recorded, all with 
CCD detectors. A double-chord occultation was observed 
in Hawaii in October 2009 and involved 2002 TX300 and 
aR = 13.4 UCAC2 star (jElliot et al. 2009)) . This was the 
first time ever that a stellar occultation by a TNO other 
than Pluto or Charon was recorded. In February 2010, 
Varuna was observed to occult a V = 11.1 star from the 
north-east of Brazil. One chord was obtained. A solution 
for the size, albedo, and shape was derived after combin- 
ing the positive with a near negative chord and astromet- 
ric and photometric ground-based data (Ortiz et al., in 
prep.). Our group subsequently observed a multi-chord oc- 
cultation of a R = 16.9 star by Eris in November 2010 
QMaury et al. 2010| |Sicardy et al. 2011b| ), a one-chord oc- 
cultation of a R = 18.2 star by 2003 AZg4 in January 2011 
( |Braga-Ribas et al. 2011a Braga-Ribas et al., in prep.), a 
multi-chord occultation of a R = 18.4 star by Makemake 
in April 2011 (Ortiz et al., in prep.) and a multi-chord 
occultation of a R — 15.7 star by Quaoar in May 2011 
(Braga-Ribas et al., in prep.). The MIT group also reported 
at the 218th meeting of the AAS in 2011 the observation 
of a one-chord occultation by Quaoar in February of 2011 
(|Person et al. 2011ISallum et al. 201 ip . All seven occulta- 
tions were predicted in this work. Our group observed them 
all (except 2002 TX300 and the February occultation by 
Quaoar) with 40 cm to 60 cm size telescopes and also with 
the NTT and VLT in Chile, Brazil, and Uruguay. 

The main difficulty in deriving reliable predictions for 
the stellar occultation of TNOs is their small apparent di- 
ameters, even in the case of larger TNOs. Another cause is 
the lack of accurate orbital elements, resulting in ephemeris 
errors as large as a few hundreds of milli-arcseconds (mas). 
Predictions based solely on published catalog positions 
such as the USNO B1.0 (jMonet et al. 2003P or the 2MASS 
(|Cutri et al. 20 03) usually fail because of the poor preci- 
sion (or lack) of proper motions and the relatively large 



M. Assafin et al.: Candidate stellar occultations by large TNOs 



3 



zonal and/or random errors in their positions (100 mas - 
200 mas) . Even individual positions for fainter stars in the 
UCAC2 catalog (Zacharia s et al. 2004[l may need correc- 
tions as large as 70 mas. 

There are ways to overcome these ephemeris and star 
position problems. The ephemeris of TNOs must be offset 
so as to get realistic TNO positions at predicted occultation 
dates. The straightforward way to do this is to regularly ob- 
serve TNOs and derive astrometrically the ephemeris offsets 
as a function of time. As successive positive occultations 
are collected, their ephemerides can be radically improved 
down to a few mas. Even a one-chord positive occultation 
helps to improve the ephemeris, since the apparent diame- 
ters of TNOs are smaller than 30 mas. In the case of star po- 
sitions, one strategy is to select possible occultations based 
on the positions given in any arbitrary astrometric catalog, 
then perform follow-up observations to improve the star 
position and (after applying offsets to the TNO ephemeris) 
delineate the shadow path. 

However, this approach to updating star positions has 
an intrinsic disadvantage in the case of TNOs. One usu- 
ally wishes to start searching for occultations from an as- 
trometrically sound catalog in order to select only those 
truly probable events. The shorter the list of candidates, 
the greater the chance of conducting the follow-up obser- 
vational programme on dedicated telescopes in due course. 
Since the highest quiality astrometric catalogs, such as the 
UCAC2 for example, are more complete for the brighter 
magnitudes, then stars fainter than say R = 15 - 16 will 
always be left out of the candidate search. This is far from 
ideal in the case of TNOs. In contrast to Pluto occultations, 
for TNOs we obtain good contrast in the light curves using 
very modest instruments for stars as faint as R = 18.0, as 
verified in the recent successfully recorded TNO occulta- 
tions reported here. This is due to the TNO magnitudes 
being much fainter than the star brightness and because 
of the recent developments in CCD detectors. Choosing 
the USNO B1.0 or other similar catalogs does not help. 
The simple addition of fainter stars that do not have the 
minimum required catalog position precision introduces the 
problem of enlarging the follow-up list with low quality tar- 
gets. This wastes valuable telescope time in follow-up pro- 
grams and may lead to the running out of nights to com- 
plete the scheduled objects. 

A more suitable strategy for TNOs is to derive local as- 
trometric star catalogs with sufficient position precision, of 
say 50 mas at least, for the time span of the occultations, 
for stars in the magnitude range R = 13 - 19. In this way, 
we match the required position precision in the search, pre- 
serving faint stars without discarding bright objects. The 
addition of astrometrically trusted faint stars to the follow- 
up list, avoiding poor quality targets in turn improves the 
chances of finding more suitable candidates for TNO occul- 
tations, owing to the increase in the star density in the sky 
path. 

The first noteworthy consistent studies of 
Pluto stellar occultation predictions were those of 



Mink fc Klemola (1985)| for the pe riod 1985-1990 and of 
McD onald fc Elliot (2000a 2000b)| for the years 1999- 
2009. Two common drawbacks of these works were the 
astrometric precision of only about 0'.'2 and the lack of 
stellar proper motions leading to uncertainties of the order 
of the Earth radius in the predicted shadow paths. In 



addition, these earlier predictions were degraded by the 
poorer precision of the older Pluto ephemerides. 

To overcome the problems faced by these first prediction 
efforts, we used more modern facilities, namely the most 
up-to-date CCD detectors, reference catalogs, and astro- 
metric procedures, and carried out an observational pro- 
gram at the ES02p2/WFI instrument. We derived local 
star catalogs with precise positions for determining accu- 
rate predictions of stellar occultations by Pluto and its 
satellites Charon, Nix, and Hyd ra for the period 20 08- 
2015. This work was published in Assafin et al. (2010) and 
is hereafter referred to as PI. Here, we extend that work 
to the predict the stellar occultations of the large TNOs 
Eris, Haumea, Ixion, Makemake, Orcus, Quaoar, Sedna, 
Varuna, 2002 TX 300 , and 2003 AZ 84 for 2008-2015, using 
the same observational and astrometric procedures. The 
local astrometric catalogs obtained are complete down to 
about R = 19, have positional errors of 40 mas or smaller, 
have stars with computed proper motions, and are in the 
UCAC2 frame. The catalogs cover the TNO sky paths with 
a 30' width in declination. They can be very useful for the 
astrometry of small CCD fields for the follow-up of selected 
candidate stars, and for refining the ephemeris offsets of 
TNOs for tracking the shadow path across the Earth. The 
catalogs are also suitable for deriving the photometric prop- 
erties of calibration and occultation stars. 

In Sect. 2, we describe the ES02p2/WFI program and 
observations. The astrometric treatment is given in Sect. 3. 
In Sect. 4, we present the derived catalogs of star positions 
along the TNO sky paths. In Sect. 5, we describe the deter- 
mination of ephemeris offsets and the candidate star search 
procedure. The predictions of the stellar occultations by the 
large TNOs are finally presented in Sect. 6. An overview is 
presented in Sect. 7, including discussions about the exter- 
nal comparison with actual occultation results and error 
budget estimates, and comments about astrometric follow- 
up programs for TNO ephemeris improvements. 



2. The ES02p2/WFI TNO program and 
observations 

The TNOs of the ES02p2/WFI program were selected 
according to their body size, which is usually related to 
their apparent sizes on the sky. Given the astrometric 
uncertainties in the star positions and the error in the 
TNO ephemeris, larger apparent diameters correspond to a 
higher probability of a successful prediction. Another selec- 
tion criteria was the star field. We preferred to select TNOs 
currently orbiting dense star regions close to the Galactic 
plane at Galactic longitudes of interest. The TNOs of in- 
trinsically high science value were also selected, even if the 
star density along their sky path was low. In these cases, 
the advantage of using the large FOV of the WFI was even 
clearer. 

Observations were made at the ESO 2p2 telescope (IAU 
code 809) using the Wide Field Imager (WFI) CCD mosaic 
detector. We used the broad-band R filter ESO#844 with 
A c = 651.725 nm and AA = 162.184 nm (full width at half 
maximum). The exposure time was 30 s . In general, S/N ra- 
tios of about 200 were reached for objects of R = 17 without 
saturating bright (R = 13-15) stars. The limiting magni- 
tude was about R = 21, and the completeness level about 
R = 19.0. The seeing varied between 0'.'6 and l'/5, being 



M. Assafin et al.: Candidate stellar occultations by large TNOs 



Table 1. Approximate {a, 6) limits and sky path periods 
covered by ES02p2/WFI mosaics for each TNO. 



TNO 


(a,S) 


(ICRF) 


(a,6) 


(ICRF) 


Time 


span 




h 


m 


O 




h 


m 


O 




yr 


Eris 


01 


38 


-04 


48 


01 


42 


-02 


58 


2008- 


2015 


Haumea 


13 


40 


+18 


57 


14 


01 


+17 


32 


2009- 


-2015 


Ixion 


16 


55 


-23 


58 


17 


28 


-27 


15 


2009- 


-2015 


Makemake 


12 


30 


+28 


25 


12 


52 


+25 


37 


2009- 


-2015 


Orcus 


09 


40 


-05 


09 


10 


04 


-08 


42 


2008- 


-2015 


Quaoar 


17 


10 


-15 


25 


17 


48 


-15 


34 


2008- 


-2015 


Sedna 


03 


27 


+06 


28 


03 


41 


+07 


17 


2009- 


-2015 


Varuna 


07 


30 


+25 


45 


08 


03 


+26 


35 


2009- 


-2015 


2002 TX300 


00 


32 


+26 


59 


01 


07 


+32 


43 


2008- 


-2015 


2003 AZ S4 


07 


36 


+ 12 


17 


08 


02 


+09 


58 


2009- 


-2015 



Approximate geocentric (a, S) limits and corresponding sky 
paths periods covered for each TNO. Paths 2008-2010 were 
observed in the September and October 2007 runs, and 
2009-2015 ones were observed in the 2009 runs at ES02p2. 
Fig. [T] illustrates the actual sky path for Varuna. 



typically 1". The telescope was shifted between exposures 
in such a way that each star in the path of the TNO was 
observed at least twice in different CCDs. 

The program was completed in seven runs, in April, 
June, September, and October 2007, May/June 2008, and 
February and December 2009. Observations spanned each 
TNO sky path from 2008 to 2015; in six cases, this path 
started in 2009, since their WFI observations only took 
place in 2009 and telescope time was limited. Table[T]shows 
the approximate limits of the sky paths observed for each 
TNO, hence indicates the period covered. Fig. [T] illustrates 
the sky path covered by Varuna for 2009-2015. A total of 
1167 WFI mosaics or 9336 individual CCD frames were 
acquired for science, resulting in about 316 GB of photo- 
metrically calibrated processed images. The two WFI field 
distortion patterns derived in PI were used here. The pat- 
tern associated with the September 2007 run was used in all 
earlier runs; the October 2007 pattern was used in all sub- 
sequent runs. The used instruments and observation proce- 
dures arc described in detail in Sect. 3 of PI. 



3. Astrometry 

We briefly describe the reduction procedures. All CCD 
images were corrected by overscan, ze romean, flatfield, 
and bad pixels with IRAF ( |Tody 1993]) using the esowfi 
(|Jones fc Valdes 2000|) and mscred (jValdes 199 8) packages. 
Using the PRAIA package (Assafin 2006), the astrometric 
treatment for each TNO consisted of three distinct steps. 

An assigned field distortion pattern was first applied to 
each CCD in the WFI mosaics for each run (see details in 
Sect. 4.1 of PI). 

In the second step, astrometry was performed over 
the individual CCDs, where the (x, y) measurements 
were corrected by the pre-determined field distortions. 
Magnitudes were obtained from PSF photometry and cal- 
ibrated with respect to the UCAC2. The (x, y) mea- 
surements were performed with two-dimensional circularly 
symmetric Gaussian fits. The estimated (x, y) measurement 
errors from Gaussian fits were about 5-10 mas for R < 16, 
then increased to 18 mas for R < 18, to 20 mas for 18 < R 
< 20, and to 30 mas at fainter magnitudes. Fig.[5]shows the 
distribution of (x, y) errors as a function of R magnitude 




7.7 7.8 7.S 

RA(hours) 



Fig. 1. Sky region covered by the ES02p2/WFI CCD mo- 
saic observations for Varuna. Years 2009-2015 follow from 
bottom to top. The continuous line is the geocentric sky 
path of Varuna. Each dashed form represents the 30' x 30' 
area covered by one single WFI mosaic. 



- a 

- $ 



*- CD 

E 



Fig. 2. (x, y) measurement errors as a function of R mag- 
nitude from all treated CCDs. Values are averages in 0.5 
magnitude bins. 



for all measurements of all TNOs. The UCAC2 was the 
prime reference catalog, except for Eris. In this case, the 
2MASS catalogue was used as the reference frame owing to 
the insufficient number of UCAC2 stars in the FOV. This 
is not an important change, as in the last step, all frames 
were consistently corrected to the UCAC2 frame. In steps 1 
and 2, we used the six constant polynomial model to relate 
the (x, y) measurements to the (X, Y) tangent plane co- 
ordinates. The position mean errors in the (a,S) individual 
CCD solutions are listed in Table [2j the average number of 
reference stars and number of CCD frames are also given 
for each TNO set of observations. Detailed information on 
the procedures in this step are found in Sect. 4.2 of PI. 

In the final step, all positions coming from the differ- 
ent CCDs and mosaics were combined to produce a global 
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Table 2. Astrometry of individual CCD frames of WFI 
mosaics. 



Table 4. Multiplicity flags for WFI global mosaic positions. 





Mean errors 


No. of 


No. ref. 


TNOs 


a(AacosS) 


a(AS) 


frames 


stars/frame 




mas 


mas 






Eris 


97 


95 


576 


27 


Haumea 


51 


50 


784 


8 


Ixion 


54 


53 


976 


108 


Makemake 


48 


51 


800 


7 


Orcus 


58 


55 


1528 


14 


Quaoar 


57 


57 


904 


121 


Sedna 


54 


57 


576 


11 


Varuna 


49 


53 


848 


31 


2002 TX300 


47 


50 


1184 


20 


2003 AZ 84 


50 


50 


1160 


48 



TNOs 


fO 


fl 


f2 


f3 


f4 


f5 




% 


% 


% 


% 


% 


% 



Mean errors are the standard deviations in the (O — C) resid- 
uals from individual (a,8) solutions with the UCAC2 cata- 
log, except for Eris in which case the 2MASS catalogue was 
used as the reference frame (see text for details). The six 
constant polynomial model was used to relate the (x, y) mea- 
surements to the (X, Y) tangent plane coordinates. The (x, 
y) centers were pre-corrected by the field distortion patterns 
(see details of these pre-corrections in Sect. 4.1 of PI). 

Table 3. Global astrometric solution for WFI CCD mo- 
saics. 







(q,<5)-UCAC2 




UCAC2 




before G.S. 


after 


G.S. 


Star 


TNOs 


Oct 




Oct 


C(5 


No. 




mas 


mas 


mas 


mas 




Eris 


62 


63 


49 


54 


232 


Haumea 


45 


44 


46 


46 


617 


Ixion 


66 


65 


53 


52 


6229 


Makemake 


41 


43 


46 


49 


523 


Orcus 


80 


73 


56 


54 


842 


Quaoar 


96 


91 


56 


55 


5996 


Sedna 


50 


53 


53 


53 


396 


Varuna 


53 


57 


49 


52 


1793 


2002 TX300 


58 


64 


46 


50 


1617 


2003 AZ 84 


56 


56 


49 


49 


3706 



The (a a , as) refer to the observed minus UCAC2 positions 
computed before and after the global solution (G.S.) proce- 
dure. Last column is the number of used UCAC2 reference 
stars. All values are averages over the results for each sky 
path year. For more details about the mosaic global solution, 
see Sect. 4.3 of PI. 



solution for each sky path year, and final {a,S) star posi- 
tions were obtained in the UCAC2 system. In this step, all 
positions were projected onto the tangent plane and a com- 
plete polynomial model of the third degree was used to fix 
all the coordinates in the UCAC2 frame. Table |3] lists the 
mean errors of UCAC2 stars before and after the WFI mo- 
saic global solutions. The listed values are averages taken 
over all sky path years. In addition to eliminating zero- 
point errors present in the positions in the individual CCD 
frames, the overall standard deviations with respect to the 
UCAC2 catalogue were generally smaller, as can be seen in 
Table [3] Detailed information about the global reduction 
procedures are given in Sect. 4.3 of PI. 

In this final step, multiple entries coming from distinct 
CCD frames for the same objects are resolved (see Sect. 4.4 
of PI for details). Table 0] displays the statistics for resolved 
multiple entries. The listed values are taken over the global 



Eris 


86.1 


4.3 


1.3 


2.8 


5.4 


0.1 


Haumea 


93.1 


0.9 


0.4 


0.7 


4.8 


0.1 


Ixion 


93.8 


0.3 


1.3 


2.6 


2.0 


0.0 


Makemake 


90.3 


1.9 


0.6 


1.2 


5.9 


0.1 


Orcus 


82.3 


9.1 


3.2 


1.9 


3.5 


0.0 


Quaoar 


92.9 


1.8 


1.0 


1.7 


2.6 


0.0 


Sedna 


92.3 


3.0 


0.9 


1.0 


2.8 


0.0 


Varuna 


93.9 


1.0 


0.5 


1.2 


3.4 


0.0 


2002 TX300 


94.2 


2.0 


1.0 


0.4 


2.3 


0.1 


2003 AZ 84 


95.6 


0.3 


0.2 


1.1 


2.7 


0.1 


All TNOs 


93.5 


0.8 


1.1 


2.2 


2.3 


0.1 



Multiplicity flag statistics for each and all TNOs. 
Percentages refer to the total number of mosaic positions 
for each TNO and for all TNOs together (last line). No flag 
(good astrometry) means no multiple entries within l'.'5 for 
the same object after mosaic solutions. Flag cases fl and 
f2 apply only to UCAC2 and 2MASS stars, when only en- 
tries assigned to a catalog were used. Flag fl means that 
more than one entry was used (more than one position en- 
try was averaged), whereas flag f2 indicates that only one 
entry was used. Flags f3, f4, and f5 apply only to field stars 
and mean that only one entry was selected according to one 
of three criteria, in order of priority: a) highest number of 
used common individual CCD positions (flag f3); b) least 
(x, y) measurement error (flag f4); c) brightest R magnitude 
(flag f5). 



solutions from all sky path years. Each final object position 
is assigned a flag that indicates its multiplicity status (see 
Table 2] for details). In all, more than 93 % of the final 
positions come with no multiple entries (flag = 0). 

Proper motions were extracted directly from the 
UCAC2 catalog. For the other stars, proper motions were 
computed for each star using the 2MASS and USNO B1.0 
catalogs as a first epoch. Although these two catalogs may 
have random/systematic errors at the 100-200 mas level, 
they are still the most suitable choice for this purpose. The 
proportion of UCAC2, 2MASS-based, and USNOBl-based 
proper motions of stars are given in Table [3] Details about 
the computations can be found in Sect. 4.5 of PI. The high- 
est percentage of stars without proper motions was verified 
for Ixion, namely 62 %, which is twice the average frac- 
tion (about 30%). Ixion is now crossing the plane of the 
Milky Way (see Table [TJ. One possible explanation is that 
the WFI magnitude limit is fainter than that of the USNO 
B1.0 (see star distribution according to magnitude in Fig. [3] 
next in the text) and many faint stars are present in this 
region of the sky. 

A complete description of the astrometric and photo- 
metric procedures used here can be found in Sect. 4 of PI, 
where data for the same WFI CCD mosaic underwent rig- 
orously the same treatment. 



4. The catalogs of star positions of each TNO for 
2008-2015 sky paths 

For each TNO, we produced star catalogs following their 
respective sky paths up to 2015. The paths begin in 2008 
for Eris, Orcus, Quaoar, and 2002 TX300, and in 2009 for 
the other TNOs. The star catalogs consist of the mean (a, 6) 
positions in the ICRS (J2000), proper motions, R magni- 
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Table 5. Proper motion computations from 2MASS, 
USNO B1.0, and ES02p2/WFI global mosaic star posi- 
tions. 



TNOs 



UCAC2 
p.m. 



2MASS 
p.m. 



USNO B1.0 
p.m. 



No 
p.m. 





% 


% 


% 


% 


Eris 


8.4 


35.5 


32.9 


23.2 


Haumea 


6.3 


24.2 


41.6 


27.9 


Ixion 


1.7 


22.3 


13.5 


62.5 


Makemake 


6.4 


25.2 


38.8 


29.6 


Orcus 


15.6 


50.3 


22.1 


12.0 


Quaoar 


5.9 


41.3 


21.6 


31.2 


Sedna 


12.2 


49.0 


22.8 


16.0 


Varuna 


10.7 


34.3 


34.7 


20.3 


2002 TX 300 


18.2 


48.0 


24.6 


09.2 


2003 AZ 84 


10.9 


27.8 


34.5 


26.8 



UCAC2 stars with proper motions (directly extracted from 
UCAC2), stars with proper motions based on the 2MASS 
and USNO B1.0, and stars for which no proper motion could 
be computed. Percentages are relative to the total number 
of stars for each TNO. 



Table 6. Star catalogs for each TNO sky path for 
2008/2009-2015. 



TNO 


No. stars 


a a 


<rs 


R magnitude 


R magnitude 






mas 


mas 


completeness 


limit 


Eris 


25,685 


33 


37 


19.0 


21.0 


Haumea 


77,860 


41 


36 


19.5 


21.5 


Ixion 


3,452,835 


41 


29 


19.5 


21.5 


Makemake 


64,487 


42 


39 


19.5 


21.5 


Orcus 


56,042 


40 


39 


18.5 


20.5 


Quaoar 


1,150,375 


41 


37 


18.5 


20.5 


Sedna 


27,464 


33 


31 


19.0 


21.0 


Varuna 


140,041 


38 


39 


19.0 


21.0 


2002 TX 300 


82,746 


33 


37 


18.0 


20.0 


2003 AZ 84 


278,720 


33 


34 


19.0 


21.5 



Number of catalog stars per TNO, position error for faint 
stars with R at completeness magnitude (errors are esti- 
mated from the standard deviation of contributing individ- 
ual CCD positions about the final catalog star positions), R 
bandpass magnitude completeness, and magnitude limit. 



tudes (also J, H, and K in the case of 2MASS stars), mean 
epoch of observations, position error at the mean epoch 
of observation, and the magnitude error estimates. In all, 
5,356,255 stars were catalogued in the UCAC2 frame. The 
mean epoch is approximately 2007.75 for the catalogs of 
TNOs beginning in 2008, and about 2009.5 for the others. 
Magnitude completeness is about R = 19. The magnitude 
limit is about R = 21. Position error is smaller than 40 mas 
for stars up to magnitude R = 19, and 30 mas up to R 
= 17. The catalogs are available in electronic form at the 
CDS. 

The catalogs are divided in terms of year. Stars that 
had multiple entries within l'/5 in the global mosaic so- 
lutions are nagged (see Table H]) . The R magnitudes from 
PSF photometry were calibrated in the UCAC2 system, 
thus magnitude zero-point errors up to 0.3 are expected 
for R > 17. The position error is estimated from repeated 
observations, by evaluating the standard deviation (mean 
error) in the contributing individual CCD positions about 
the final catalog star positions (last step of global mosaic 
solution - see Sect. 3). By default, multiple-entry nagged 
stars have no position error estimates. Infrared magnitudes 
(and errors) were extracted from the 2MASS catalog. Error 
estimates for R magnitudes come from the standard devi- 
ation about the mean for individual CCD frames. The sky 
coverages of the catalogs have widths of at least 30' in dec- 
lination and lenghts of a few degrees along the geocentric 
sky path of each TNO. The proportion of UCAC2, 2MASS, 
and USNO B1.0 stars in the WFI catalogs may be verified 
from Table [5] 

For each TNO, Table [5] lists the total number of catalog 
stars per year, the R bandpass magnitude completeness, 
magnitude limit, and average position errors for faint stars 
with R at the magnitude completness. Fig. [3] shows the star 
distribution per R magnitude for all catalogs. Fig. [4] plots 
the position error of all catalog stars as a function of R 
magnitude. Values were averaged over 0.5 magnitude bins. 




(%) sjeis -on 

Fig. 3. Star distribution per R magnitude for all catalog 
stars. It illustrates the R magnitude completeness and limit 
of the catalogs. Percentages were computed over 0.5 mag- 
nitude bins for all stars in all catalogs. 



5. Ephemeris offsets and search for candidate 
stellar occultations 

To achieve the highest accuracy in predicting stellar oc- 
cultations, it is necessary to pin down the actual paths 
of the TNOs in the sky. This can be done by correcting 
their orbits by offseting their ephemerides in both right 
ascension and declination. We derived average ephemeris 
offsets for all TNOs in the ES02p2/WFI program. Table [7] 
displays the ephemeris offsets information for each TNO. 
The TNO observed positions used to derive the offsets 
were in the UCAC2 system. They came from a parallel 
astrometric program of CCD observations carried out dur- 
ing 2007-2010 at a variety of 0.6 m to 2.2 m telescopes 
in Chile (ES02p2/WFI), Brazil (Laboratorio Nacional de 
Astrofi'sica - LNA) , and France (Pic du Midi; Observatoirc 
Haute de Provence). The obtained offsets are provisional 
in the sense that the aforementioned program is still ongo- 



M. Assafin et al.: Candidate stellar occultations by large TNOs 



7 



□ right ascension 
x declination 



n □ □ □ □ □ L 



10 12 14 



20 22 24 



R magnitude 



Fig. 4. Position mean errors as a function of R magnitude 
for all catalog stars. Position errors are estimated from the 
standard deviation in the contributing individual CCD po- 
sitions about the final catalog star positions (last step of 
global mosaic solution - see Sect. 3). Values were computed 
over 0.5 magnitude bins. 

Table 7. TNOs offsets with respect to JPL ephemerides. 



Using the astrometric catalog of stars in the path of 
the sky through which the TNOs will pass up to 2015, 
we then searched for stars that could be occulted by the 
respective TNOs. To achieve this, the whole star catalogs 
were cross-correlated with the JPL ephemerides of the bod- 
ies generated each minute for the whole period between 
2008/2009 to 2015. The JPL ephemerides were offset prior 
to the search using the quantities given in Table [7] If the 
apparent distance in the sky plane between the star and 
the body is less than an Earth apparent radius plus 50 mas 
of safe margin, then a potential occultation was identified 
and all astrometric and geometric data relevant to the pos- 
sible event were computed and stored. In addition to as- 
trometric and photometric data, for each candidate star, 
we computed and stored the minimum apparent geocentric 
distance d, the central instant of closest approach to, the 
shadow velocity v across the Earth, the position angle P/A 
of the shadow path, and the local solar time LST at the 
sub-planet point. These geometric quantities are defined 
and calculated as explained in Sect. 7 of PI. 



6. Predictions of stellar occultations by ten large 
TNOs 

Following the procedure described in Sect. 5, we iden- 
tified potential candidate stars for occultations by Eris, 
Haumea, Ixion, Makemake, Orcus, Quaoar, Sedna, Varuna, 
2002 TX300, and 2003 AZg 4 . The adopted search radius was 
0'.'335 - i.e. 285 mas (approximately the apparent Earth ra- 



Observed - JPL ephemeris 



No. 



JPL dius at 31 AU, i.e., at Pluto for 2008-2015) plus a safe 





Aacos(5 


AS 




E s 


Nights 


epheme 




mas 


mas 


mas 


mas 




versio 


Eris 


-020 


-068 


084 


055 


09 


28 


Haumea 


-090 


-098 


060 


081 


21 


37 


Ixion 


+020 


+191 


092 


093 


15 


09 


Makemake 


-039 


-165 


060 


058 


20 


40 


Orcus 


-087 


+001 


059 


046 


19 


13 


Quaoar 


-144 


-078 


098 


043 


22 


17 


Sedna 


-043 


-160 


040 


080 


05 


17 


Varuna 


-045 


+005 


015 


015 




21 


2002 TX300 


-084 


-053 


040 


040 




14 


2003 AZ 84 


+232 


+006 


003 


015 
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Average "observed minus JPL ephemeris" position offsets. 
Errors come from the standard deviation about the mean of 
individual nightly averages. TNO positions in the UCAC2 
system come from a parallel astrometric program of CCD 
observations carried out during 2007-2010 at a variety of 
0.6 m to 2.2 m telescopes in Chile, Brazil, and France. Each 
TNO ephemeris version is displayed. They were extracted 
using the JPL Horizons service (Giorgini et al. 1996| ). The 
listed offsets are provisional in the sense that the parallel 
program is still ongoing and further updates are expected. 
For three TNOs (Varuna, 2002 TX300, 2003 AZ 84 ), offsets 
came from fits to actual observed stellar occultations (see 
discussion in Sect. 7). 



ing and further refinements and updates are expected. For 
three TNOs (Varuna, 2002 TX 300 , 2003 AZ 84 ), we derived 
the offsets from fits to actual observed stellar occultations 
(see discussion in next Sect. 7). The offset error estimates 
come from the standard deviation about the mean of indi- 
vidual nightly averages. They do not depend upon the par- 
ticular ephemeris errors themselves, which were published 
by the JPL. 



even so, be visible right above the horizon from places still 
in the dark near the Earth terminator. For each TNO, all 
relevant astrometric, photometric, and geometric informa- 
tion for each potential event found is available in electronic 
fornQ. Table |H] lists a sample of predictions for 2002 TX300, 
Varuna, Eris, 2003 AZg 4 , Makemake, and Quaoar. These 
very predictions were used to prepare campaigns actu- 
ally resulting in successful recordings of stellar occultations 
by our group. The exception was 2002 TX300, for which 
no observations could be done at Hawaii (the predicted 
path), because all telescopes were occupied, must of them 
having been scheduled to observe the impact of NASA's 
Lunar Crater Observation and Sensing Satellite (LCROSS; 
IGarner fc Dun bar 2009), which occurred about one hour 
after the predicted central instant of occultation. Prediction 
tables contain the date and instant of stellar occultation 
(UTC), both the ICRS (J2000) star and TNO estimated 
coordinates at the event date, the closest apparent geocen- 
tric distance between the star and TNO, the position angle 
of the shadow across the Earth (clockwise, zero at North), 
the velocity in km s , the distance to the Earth (AU), lon- 
gitude of the sub-solar point, local solar time, applied JPL 
ephemerides offsets in (RA, DEC) for the central instant, 
catalog proper motion and multiplicity flags, estimated star 
catalog position errors, proper motions, and magnitudes 
R*, J*, H*, and K*. Magnitudes were normalized to a ref- 
erence shadow velocity of 20 km s^ 1 by 



M* = M + 2.5 log 



10 



20 km s 



anonymous ftp to cdsarc.u-strasbg.fr 
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Table 9. Predictions of stellar occultations by ten large 
TNOs for 2008-2015. 



close-by faint objects, and other elements in the images may 
prevent problems in the recording and photometry of oc- 
cultation observations. When it comes to an event, charts 



TNO 


2008 


2009 


2010 


2011 


2012 


2013 


2014 


Eris 


2 


1 


1 








1 





Haumea 


X 


3 


4 


6 


4 


7 


7 


Ixion 


X 


102 


193 


288 


352 


425 


309 


Makemake 


X 


4 


3 


5 


3 


1 


4 


Orcus 


4 


3 


6 


1 


3 


1 


2 


Quaoar 


51 


26 


27 


90 


176 


53 


67 


Sedna 


X 


1 





2 


1 


2 


3 


Varuna 


X 


9 


11 


11 


8 


14 


7 


2002 TX300 


1 


4 


3 


3 


2 


5 


1 


2003 AZ S4 


X 


18 


18 


19 


16 


9 


15 



Number of predicted events per year for each TNO. 

The reference velocity of 20 km s _1 is typical of events 
around Pluto opposition. In slow events, M* may become 
one to two magnitudes brighter than the actual star mag- 
nitude, indicating that the observation is feasible after all, 
because we can use longer integration times, and conse- 
quently reach reasonably high S/N, without any losses of 
spatial resolution in diameter measurements. For the same 
reasoning, we may also probe the atmosphere altitudes in 
the light curves, despite the faintnesses of the stars. As 
mentioned in the introduction, in the case of TNOs, oc- 
cultations involving stars as faint as R = 18.0 have been 
successfully recorded with 50 cm telescopes. 

Table [9] displays the total number of predicted events 
for each TNO for the period 2008/2009-2015. 

7. Discussion 

We have presented predictions for stellar occultations by 
the large TNOs Eris, Haumea, Ixion, Makemake, Orcus, 
Quaoar, Sedna, Varuna, 2002 TX 300 , and 2003 AZ 84 for 
2008/2009-2015. We applied the same observational and 
astrometric procedures as a pre vious study dedicated to the 
Pluto system and published in Assafin et al. (2010)| Many 
issues were addressed in that work, such as body ephemeris 
offset determination, catalog zero-point position errors and 
field-of-view size, long-term predictions and stellar proper 
motions, and faint- visual versus bright infrared stars, which 
also hold true here. 

The local astrometric catalogs that we obtain are com- 
plete to about R — 19 and in the UCAC2 frame. Up to R 
magnitude completeness brightness, stars have positional 
errors of 40 mas or smaller and computed proper motions. 
The catalogs cover the TNO sky paths with 30' width in 
declination and can be very useful in astrometric follow-up 
programs on telescopes/CCDs with FOVs of any size. The 
catalogs are also sources of photometric information for cal- 
ibration and occultation stars. The astrometry was made 
by the use of the PRAIA package (| Assafin 2006[) . About 
316 GB of science-alone WFI images were treated, not to 
mention calibration observations to derive field distorion 
patterns, which improved the astrometric results of WFI 
data to the 40 mas level (see Sect. 9 in PI for a detailed 
discussion about the field distortion pattern improvements 
for WFI astrometry). 

As an important by-product, the high- resolution CCD 
frames obtained (1 pixel = 0'.'238) form an image bank use- 
ful for the visual inspection of the fields of view. For in- 
stance, the analysis of small details in the calibration stars, 



exffi&ftcd from this image bank are usually posted in our 
wcbOpage, but any particular field around some desired po- 
sition can be obtained in FITS format on request to the 
au tifrs. Public access to the original images from the ob- 
served WFI CCD mosaics will soon be available through 
the^ESO Image Archive Portal. 

]§o threshold in R magnitude was established in the 
seargh for occultations, as faint R objects may turn out to 
be bright IR stars, suitable for observation at adequately 
equipped IR ground-based f acilities or by the SOFIA ob- 



servatory (|Gehrz et al. 2009]) ; H, J, and K magnitudes are 
readily available in the catalog if the star belongs to the 
2MASS. In addition, events that occur at slow shadow 
speeds of less than 20 km s _1 may become observable. No 
constraints on a geographic place were applied, as in prin- 
ciple SOFIA observations can be done from any sub-solar 
point on Earth. Events in daylight at sub-planet point were 
not excluded either, as they remained observable in the 
dark, right above the horizon, from places near the Earth 
terminator. Throughout this paper, we do not distinguish 
between past and future predictions, publishing all found 
occultations for the sky path covered (or to be covered) 
by the TNOs over the years 2008/2009-2015. We assume 
that the WFI predictions of past occultations may be useful 
as either reference for ongoing fittings of recently obtained 
light curves, an aid to deriving ephemeris offsets by compar- 
ing expected and observed central instants and C/A values, 
or an external check for the accuracy and precision of our 
WFI predictions. 

The information in the prediction tables allows us to de- 
termine the central instant and shadow path over the Earth. 
For ground-based observations, the selection of events usu- 
ally favors brighter stars and shadow paths that cross eas- 
ily accessible places for the use of mobile telescopes. Faint 
stars may be selected in the case of observation at sites 
with larger observatories/telescopes. If the star is reddish, 
observations can be made in daylight with IR detectors. 
For TNOs with a low rate of occultations per year, such 
as Eris, opportunities are very rare, so broader nets of ob- 
servers may be formed over larger areas, even if the event is 
not as promising in terms of star brightness or location. We 
leave the choice of potential events among the many pre- 
dicted occultations listed in the tables to the judgement of 
the reader. Our group publishes at our web page potential 
occultation campaigns on a yearly basis. 

The errors in our TNO sky path predictions are gov- 
erned by the internal offset uncertainties listed in Table [3 
These uncertainties are based on optical observations made 
for a period of some years. The budget ephemeris error es- 
timates given by JPL - which increase with each passing 
year - only express the random error features of the JPL 
orbit solutions, rather than represent a possible increase in 
systematic trends. As time goes by and more TNO obser- 
vations and occultations are accumulated, more accurate 
JPL ephemeris solutions are being published, although in 
the new versions no significant shifts in the TNO orbital 
sky paths are predicted for the next few years. For now, 
the adoption of a constant ephemeris offset in the predic- 
tions has the only side effect of absorbing some remaining 
systematic errors into the offset uncertainties. These un- 
certainties, on the other hand, do not increase with time 
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year m d h m s 


RA (ICRS) Dec 


C/A P/A 


v D 


R* J* H* K* 


A LST 


Ae a Aeg pm ct fg 




h m s ° ' 


mas ° 


km s" 1 AU 




° h:m 


mas mas 



2009 10 09 10 29 29 00 37 13.6205 +28 22 23.019 034 158.89 -25.29 40.61 13.5 12.6 12.2 12.2 194 23:24 -84.0 -53.0 ok uc 

2010 02 19 23 07 40 07 29 22.4672 +26 07 23.191 031 188.28 -18.66 42.75 11.8 10.3 10.0 10.0 336 21:30 -45.0 +05.0 ok uc 

2010 11 06 02 08 11 01 39 09.9421 -04 21 12.119 021 347.56 -26.32 95.74 16.8 14.3 13.7 13.5 308 22:38 -20.0 -68.0 ok 2m 

2011 01 08 06 31 48 07 43 41.8220 +11 30 23.569 094 187.02 -26.00 44.32 18.5 50.0 50.0 50.0 271 00:33 +232.0 +06.0 ok fs 
2011 04 23 01 37 59 12 36 11.3938 +28 11 10.493 174 183.80 -22.09 51.50 18.4 50.0 50.0 50.0 314 22:33 -39.0 -165.0 ok fs 
2011 05 04 02 40 45 17 28 50.8021 -15 27 42.788 046 191.77 -18.28 42.35 15.7 13.5 12.8 12.7 000 02:42 -144.0 -78.0 ok uc 

In entry order, the above samples correspond to the predictions for 2002 TX300 (October 2009), Varuna (February 2010), Eris 
(November 2010), 2003 AZ 8 4 (January 2011), Makemake (April 2011), and Quaoar (May 2011), respectively. Prediction tables 
list event date and instant (UTC), the ICRS (J2000) star coordinates at occultation, the closest apparent distance between 
star and TNO (C/A), the position angle (P/A) of the shadow across the Earth (counter-clockwise, zero at north), the velocity 
in km s _1 , the distance (D) to the Earth (AU), R*, J*, H*, and K* star magnitudes normalized to a reference shadow velocity 
(v) of 20 km s _1 (50.0 means no magnitude available), longitude (A) of the sub-solar point, local solar time (LST), (Ae a , Aes) 
JPL ephemeris offsets in (a, S) for the central instant (see Sect. 5), catalog cross-identification (uc = UCAC2, 2m = 2MASS, 
fs = field star), proper motion existence and multiplicity flags (see Sect. 3), estimated star-catalog position errors (E a , Ea), 
and proper motions (fl a , fig)- Positive/negative v means, respectively, prograde/retrograde velocities, that is TNO's geocentric 
right ascension is increasing/decreasing, respectively. The complete table sets of 2008/2009-2015 predictions for Eris, Haumea, 
Ixion, Makemake, Orcus, Quaoar, Sedna, Varuna, 2002 TX300, and 2003 AZs4 are available in electronic form at the CDS. In 
the electronic version, in addition to the star positions, we provide the TNO estimated coordinates (with ephemeris offsets 
applied) at occultation instant. 



Table 10. TNO diameters. 



TNO diameter 





km 




mas 


Eris 


2326 (+/- 


-12) 


33 


Haumea 


1324 (+/- 


167) 


36 


Ixion 


0650 (+250/ 


-220) 


22 


Makemake 


1455 (+/- 


-12) 


38 


Orcus 


0850 (+/- 


-90) 


25 


Quaoar 


1170 (+/- 


-40) 


34 


Sedna 


1200-1600 


22 


Varuna 


1003 (+/- 


-09) 


32 


2002 TX300 


0143 (+/- 


-05) 


05 


2003 AZ 84 


0910 (+/- 


-60) 


28 



dictions (ES02p2/WFI-based) and actually observed oc- 
cultations, for which P/A, central instant and C/A val- 
ues were obtained from fits to the light curves of the 
Reference c hor ds. We used six out of the seven TNO stellar oc- 
ISicardy et al. 201R2 jultations known to date, for which data is available 
Lellouch ct al. 201Q EUi ot ef &{ 2 QQC;| |Sicardv et al. 201 lb| Braga-Ribas, 2011, 

Stansbcrry et al. 200» . — \~TT- TTTT — ' j w j 

Ortiz e t al 201l f ^ lv ' comm -)- Using P/As, distances and relative speeds, 



|Lim, 



we computed both the predicted minus fitted central in- 



et al. 2010, 

|Braga-Ribas et al~20 l'iB [1 ^ anc ^ C/A differences. Differences in central instants 
IBrown 20081 are expressed in seconds and in mas (in order to be compa- 
Sicardy et al. 2010f able with C/A results), and C/A differences are given in 
Elliot ct al. 2009 both mas and kilometers. We give in Table [TT] the averages 
Mueller ct al. 2010 : or the modulus of the differences. We note that the WFI 



TNO diameters in km (uncertainties in parentheses) and 
as geocentric apparent sizes in mas for 2012.5. Values were 
obtained from previous stellar occultations and by modeling 
visible and IR observations (see references). 



for the period considered and are realistic estimates for 
the overall sky path prediction uncertainties. In the pre- 
diction tables (available in eletronic form), in addition to 
the star position, we indicate the estimated TNO position 
(with ephemeris offset applied) at the predicted occultation 
instant. This should help us to evaluate the events in the 
future as updates become available at the star position and 
on the TNO ephemeris by follow-up astrometric programs 
or by the JPL. 

A tentative way to estimate the probability of success 
in recording an occultation of a particular TNO is to con- 
sider the ratio of the apparent TNO diameter in mas to its 
ephemeris offset errors displayed in Table [7] We list TNO 
diameter values in Table [10] to support reliability estimates 
for the predictions. They are given in km and as geocen- 
tric apparent sizes in mas for 2012.5. Values were obtained 
from previous stellar occultations and by modeling visible 
and IR observations (see references in Table ITU1) . 

The ES02p2/WFI predictions can be directly com- 
pared to actually observed stellar occultations in order 
to estimate the external accuracy and precision of predic- 
tions. Table [11] displays a direct comparison between pre- 



prediction for the Eris' central instant was in strong dis- 
agreement with observations. This may indicate that much 
more accurate right ascension ephemeris offsets should be 
obtained in our follow-up program for Eris. A separate set 
of averages is provided in Table [TT] without taking Eris into 
account. The TNO 2003 AZ§4 is excluded from the analy- 
sis as the fit to the actual occultation is based on the WFI 
star position, and the current WFI predictions already uses 
the offsets derived from the occultation fitting. In all, the 
differences in C/A and central instant are consistent with 
the expected errors in the WFI star positions as a func- 
tion of magnitude and ephemeris offset error (see Table [7]) . 
The somewhat larger differences found in Table [TT1 for the 
central instant may be associated with ephemeris offset un- 
certainties in right ascension (the direction of the chords are 
usually east- west). One possible cause is geocentric paral- 
lax error (see discussion about the follow-up program and 
ephemeris offset determination in the next paragraphs). In 
all, for R = 19 stars (catalog magnitude completeness) and 
40 mas errors in the WFI positions, we may assume a bulk 
error of about 80 mas for C/A, that is dominated by the 
ephemeris offsets errors of about 70 mas. For about 40 A.U., 
this implies a shadow path uncertainty over the Earth of 
the order of 2300 km. If the ephemeris offsets can be well- 
determined to a more accurate precision than 30 mas, then 
a bulk error of 50 mas in C/A can be achieved, leading 
to a precision of about 1400 km for the WFI occultation 
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Table 11. Comparison between WFI-based and actual ob- 
served occultations: central instants and C/As 



Star 



Prediction minus Occupation 



TNO 


mag (R) 


Central instant 




C/A 






s 


mas 


mas 


km 


2002 TX 300 


12.9 


-15 


-13 


-09 


-265 


Varuna 


11.0 


-20 


-12 


-63 


-1953 


Eris 


16.9 


-604 


-229 


+39 


+2708 


2003 AZ 84 


18.2 


+0 


+0 


+0 


+0 


Makemake 


18.4 


+196 


+116 


+40 


+1494 


Quaoar 


15.7 


+130 


+78 


-44 


-1351 


average (moduli) 




193 


90 


39 


1554 


average (Eris out) 




90 


55 


39 


1266 



Comparison between central instants and C/As from WFI- 
based (Table [H| and actual observations (fit to data) for 6 
(out of seven) successfully recorded TNO stellar occultations 
to date for which data is available. Using P/A, distance and 
the relative speed, central instant differences are expressed in 
seconds and in mas, and C/A differences in mas and in km. 
Averages for the modulus of the differences are computed 
with and without Eris. 2003 AZs4 is excluded from the anal- 
ysis as the fit to the actual occultation already uses the WFI 
star position, and the current WFI predictions already uses 
the offsets derived from the occultation fitting. 



path predictions, which is similar to the obtained values 
displayed in Table [TO 

Long-term predictions are the only way of spotting po- 
tential stellar occultations in terms of favourable geographic 
access, observatories and instrumental facilities, etc. Once 
an event is chosen, then a complex international campaign 
is initiated, with astrometric follow-up of the TNO and star 
being mandatory. The only reliable long-term predictions 
of stellar occultations that use solely all-sky star catalogs 
are those based on the UCAC2. This limits the search for 
candidates to about R = 16 stars. The goal of our WFI pro- 
gram was to extend the UCAC2 position precision to stars 
as faint as about R = 19. The availability of data for such 
faint stars with about 40 mas position precision improves 
the chances of finding more suitable candidates for TNO 
occultations, owing to the increase in the number of suit- 
able stars in the sky path and, consequently, the frequency 
of appulses. Owing to the recent developments of solid state 
detectors, and because TNOs are usually much fainter than 
the occulted stars, we can obtain good contrast in the light 
curves, even using very modest instruments. This has been 
verified for R = 17 - 18 stars in the occultations by Eris 
and 2003 AZ 84 ([Maury et al. 2010[ |Sicardy et al. 2011b 
Braga-Ri bas et al. 201 la| Braga-Ribas et al. 2011, in 



prep.), observed with 40 cm and 60 cm telescopes in Chile. 

The sky paths covered by the star catalogs were deter- 
mined by the chosen time span of our predictions (2008- 
2015) in a compromise between the availability of WFI tele- 
scope time and the number of selected TNOs. The upper 
limit coincides with NASA's New Horizons spacecraft ar- 
rival at Pluto in 2015. Given this long time coverage of 
about eight years, the computation of proper motions for 
the stars was also of concern. That is why we used as first 
epoch data the best available astrometric resources at the 
time of completion of this work. In all, the comparison made 
with actual events (Table [TTj) and the successful WFI pre- 
dictions specifically for the faint stars of Eris, 2003 AZ§4 
and Makemake occultations indicate that the WFI astro- 



metric solution is sactisfatory in the context of long-term 
predictions. 

Owing to the lack of observations of trans-neptunian 
"objects for ephemerides, to the few positions sometimes 
derived without the appropriate astrometric care, and to 
the relatively short orbital arc length coverage, the ac- 
tual position of TNOs can differ significantly from their 
ephemerides. To refine the orbits of large TNOs in the short 
and long runs, by regularly obtaining precise astrometric 
positions, a parallel program involving many telescopes of 
60 cm to 2 m class in Chile, France, and Brazil was started 
in 2007 (see further details in Sect. 5). Using the first re- 
sults of this program for the predictions published in this 
work, we applied average fixed offsets to the current JPL 
ephemerides of the TNOs in the list (Table HJ). In some 
cases, offsets reached almost 200 mas (about one Earth ra- 
dius at 30 AU) or more. 

Given the time span of 2008-2015, we expect that the 
fixed corrections applied to the JPL ephemerides suffice 
for pinpointing the most interesting occultations among all 
predictions. However, once an event is selected, the shadow 
path on Earth must be further refined. From the point of 
view of TNO orbits, improvements can still be made to the 
accuracy of the WFI predictions. For instance, by observ- 
ing TNOs at their opposition and conjunctions to obtain 
the observed ephemeris offsets, one can model the geocen- 
tric parallax errors and correctly extrapolate the observed 
ephemeris offsets with time for the occultation dates. For 
a body at 30 AU in the Ecliptic, these errors can amount 
to 100 mas in right ascension for an ephemeris error of 
some few tens of thousand kilometers in heliocentric dis- 



tance ( Camargo et al. 2010 1. As successful occultations are 
recorded, the ephemeris offsets are more accurately deter- 
mined. For three TNOs (see Table [7]), we could already im- 
prove the ephemeris offsets for this work using occultation 
results. As the new occultation results are consolidated, we 
will update the ephemeris offsets and release new, higher 
quality predictions. 

An advancement that is beyond the scope of WFI pre- 
dictions, for selected occultations, is an alternative proce- 
dure for improving the shadow path in the short term by 
direct fitting an orbital arc to the observed TNO positions. 
This can be done by numerically integrating orbits or, in 
principle, by using the current ephemeris as a template 
and then fitting it to the observations following a correc- 
tion model. In the case of observations when both star and 
TNO are imaged in the same field of view, this procedure 
can yield very accurate predictions, as high precision rela- 
tive positions can be attained for the TNO with respect to 
the occultation star, once proper care is taken to correct 
for differential color refraction. 

The astrometric follow-up of the selected stars and 
TNOs on a regular basis is therefore very important, re- 
gardless of the approach taken to improve the accuracy of 
the predictions. In practice, however, the week or two be- 
fore occultation date is still decisive in pinning down the 
shadow path on Earth. Updates on the TNO ephemeris 
offsets, star positions of selected occultations, updated re- 
ports, and finding charts are available on a continuous ba- 
sis by our group at our web pagfi Among other similar 
projects conducted by other active groups in the U.S.A., 
our efforts reflect a long-term collaboration with the in- 



http:/ /www. lesia.obspm.fr/perso/bruno-sicardy/ 
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ternational campaigns coordinated by the Observatoire de 
Paris for this purpose. 
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